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ABSTRACT 

The binaries PSR Jl 141-6545 and PSR B2303+46 each appear to contain a white 
dwarf which formed before the neutron star. We describe an evolutionary pathway to 
produce these two systems. In this scenario, the primary transfers its envelope onto 
the secondary which is then the more massive of the two stars, and indeed sufficiently 
massive later to produce a neutron star via a supernova. The core of the primary 
produces a massive white dwarf which enters into a common envelope with the core 
of the secondary when the latter evolves off the main sequence. During the common 
envelope phase, the white dwarf and the core of the secondary spiral together as the 
envelope is ejected. The evolutionary history of PSR J1141-6545 and PSR B2303+46 
differ after this phase. In the case of PSR J1141-6545, the secondary (now a helium 
star) evolves into contact transferring its envelope onto the white dwarf. We propose 
that the vast majority of this material is in fact ejected from the system. The remains 
of the secondary then explode as a supernova producing a neutron star. Generally 
the white dwarf and neutron star will remain bound in tight, often eccentric, systems 
resembling PSR J1141-6545. These systems will spiral in and merge on a relatively 
short timescale and may make a significant contribution to the population of gamma 
ray burst progenitors. In PSR B2303+46, the helium-star secondary and white dwarf 
never come into contact. Rather the helium star loses its envelope via a wind, which 
increases the binary separation slightly. Only a small fraction of such systems will 
remain bound when the neutron star is formed (as the systems are wider). Those 
systems which are broken up will produce a population of high-velocity white dwarfs 
and neutron stars. 

Key words: accretion, accretion discs — binaries: close stars: evolution — stars: 
stars. 



1 INTRODUCTION 

Observations suggest that in the two binaries PSRs 
B2303+46 and Jl 141-6545, the companions of the observed 
pulsars are white dwarfs which appear to have been formed 
before the neutron stars (van Kerkwijk & Kulkarni 1999; 
Kaspi et al 2000; Manchester et al 2000). 

If the white dwarf were the remains of a former donor 
star of the neutron star, we would expect the white dwarf to 
be in a circular orbit. However B2023+46 has an eccentricity 
of 0.658, so we must conclude that the white dwarf was made 
first and the current observed eccentricity is derived from the 
kick received when the neutron star was formed. 

J1141-6545 is a non-recycled pulsar with a massive 
companion in an eccentric orbit (Kaspi et al 2000; Manch- 
ester et al 2000). If the pulsar had formed first, we would 
expect it to have been spun up (ie recycled) when the sec- 
ondary evolved into contact and transferred material onto 



it, ultimately leaving either a second neutron star or a white 
dwarf. If the companion produced a white dwarf after the 
birth of the neutron star, we would also expect the orbit 
today to be circular, whereas it has a measured eccentric- 
ity of 0.172. Hence it seems likely that in J1141-6545 the 
companion is also probably a massive white dwarf which 
was formed before the neutron star. It should be noted that 
although many of the basic properties of J1141-6545 and 
B2023+46 are similar, their orbital periods differ by a fac- 
tor of 60. Accounting for this fact is one of the problems we 
are going to address in this paper. 

We consider an evolutionary scenario for the produc- 
tion of such white-dwarf neutron-star (WD-NS) binaries. In 
this scenario, the original primary transfers its envelope to 
the secondary conservatively leaving only the helium-star 
core. This helium star then evolves, filling its Roche lobe, 
leading to a second phase of mass transfer where the en- 
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velope of the helium star is transferred to the secondary 
(which is still on the main sequence) . In the process the pri- 
mary becomes a white dwarf. In sufficiently wide binaries, 
when the secondary evolves into contact, the ensuing phase 
of mass transfer will produce a common-envelope system in 
which the white dwarf and the helium-star core of the sec- 
ondary will spiral together as the common envelope of gas is 
ejected from the system. The post-common-envelope system 
will consist of the white dwarf and the helium star in a tight 
binary. Providing sufficient mass has been transferred to the 
secondary during earlier phases in the binary evolution, the 
helium star may be massive enough to produce a neutron 
star via a supernova explosion. 

This scenario has been considered elsewhere (see for ex- 
ample Portegies-Zwart & Yungelson 1999; Brown et al 2000; 
Tauris & Sennels 2000). Here we extend the earlier work by 
considering the role of the helium star radius in the evolution 
of the system after the common envelope phase which will 
depend on the separation of the binary and the maximum 
size of the pre-supernova helium star. 

Systems similar to B2023+46 will be produced in rela- 
tively wide systems, where the helium star avoids filling its 
Roche lobe as it expands. Instead it loses mass via a wind. 
For closer binaries (ie those producing systems similar to 
Jl 141-6545), the helium star in the post-common-envelope 
system will fill its Roche lobe. In this paper, we propose 
that a phase of mass transfer follows where the helium star 
envelope is transferred to the white dwarf at highly super- 
Eddington rates. Rather than being accreted, we propose 
that this material is ejected from the system. 

In Section 2, we begin by describing a simplified version 
of the scenario described above where we neglect the radius 
of the final helium star. In Section 3 we compute the popu- 
lation of binaries produced using this simplified scenario. In 
Section 4 we consider the role of the helium-star radius on 
the binary evolution. In Section 5 we consider the inferred 
birthrate of B2023+46-like and J1141-6545-like systems. 
We also consider the timescale for the merger of the tighter 
systems and the importance of the kick velocities received by 
the binaries on formation of the neutron star. We summarise 
our results in Section 6. 



2 EVOLUTIONARY SCENARIO 

The key feature of both J1141-6545 and B2303+46 is that 
the neutron star is younger than the white dwarf. This im- 
plies that the star which produced the neutron star in each 
system must have originally been the secondary. Yet the 
primary was only massive enough to produce a white dwarf. 
Mass transfer from the primary to the secondary must have 
increased its mass sufficiently to yield a supernova. These 
requirements greatly restrict the binary parameter space 
which may produce such neutron star-white dwarf (WD- 
NS) systems. The evolutionary sequence envisaged to pro- 
duce WD-NS binaries is outlined below: 

We consider a binary of initial separation, ai, and total 
mass contained in the system Mi = Mi,i + M 2 ,i. 

The first phase of mass transfer occurs when the pri- 
mary fills its Roche lobe just off the main sequence. Dur- 
ing this phase of conservative, radiative mass transfer (early 
Case B), the primary transfers its entire envelope to the sec- 



ondary. Only the helium core of the primary remains. The 
masses of the two stars, Mi,b and M 2 ,b, after this phase of 
mass transfer are given by 



Mi, B = M He (Mi,i) = aMli 



M 2 ,b = Mi - Mi, : 



(1) 
(2) 



where a = 0.125 and b = 1.4 (van den Heuvel 1994). The ra- 
tio of the separation after this phase to the initial separation 
is given by 



fUe = — 



a B Mi-Mi 



(3) 



This helium star then evolves filling its Roche lobe lead- 
ing to the second phase of mass transfer (Case BB) where 
the primary transfers the envelope of the helium star to the 
secondary (which is still on the main sequence). In a calcu- 
lation of Case BB evolution by Delgado & Thomas (1981), 
a 2 M Q helium star initiates mass transfer when its radius 
is about 20 R0 and terminates mass transfer at a radius of 
about 50 Rq yet it is stripped down to the CO core. We 
therefore assume here that after this phase of mass trans- 
fer the primary becomes a CO or ONeMg white dwarf. The 
masses of the two stars, M^bb and M 2j bb, after this phase 
of mass transfer are given by 



Mi, BB = M CO (Mi,b) = Mwd 



M 2 ,BB =M\~ Ml,BB 



(4) 

(5) 



The ratio of the separation after the Case BB mass 
transfer to the initial separation is given by 



JU B B = 



«BB 

OSi 



MljMjj 



(6) 



The secondary will evolve into contact on its nuclear 
evolutionary timescale. If the system is sufficiently wide, 
this third phase of mass transfer will be convective Case B 
which will lead to a common-envelope phase during which 
the white dwarf and the helium-star core of the secondary 
spiral together as the common envelope of gas is ejected from 
the system. Alternatively, the a common-envelope phase will 
be initiated from radiative, Case B mass transfer, followed 
by a delayed dynamical instability (DDI) if the mass ratio is 
large enough (Hjellming 1989, Kalogera & Webbink 1996). 
The masses of the two stars after this phase, Ml,ce and 
M 2j ce, are given by 



Mi,cb = Mi, BB = M W d 



M 2 ,CE = M H e(M 2 ,BB) = OM^BB 



(7) 
(8) 



The inspiral during the common envelope phase may be 
computed by equating the change in orbital energy of the 
two stars to the binding energy of the envelope up to an 
efficiency qce, ie E cnv = oceA_E. Here E env can be written 
in the following form (Webbink 1984) 



E — 



GM 2 ,bb(M 2 ,bb - M 2 ,ce) 
Ace/ 2 {q BB )ci bb 



(9) 



where the radius of star 2 when it fills its Roche lobe is given 
by -R 2 ,bb = / 2 (<7BB)aBB, where q BB = Mi i bb/M 2 ,bb and 
f 2 (q) = 0.49g 2/3 /(0.6g 2/3 + ln(l + g 1/3 )) (Egglcton 1983). 
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Combining equation (9) with the expression for AE, after 
some rearrangement we arrive at the following expression 
for the inspiral during the common envelope phase 



Jbb^ce 



acE - 2M2,bb(M2,BB — M2,ce) _|_ 

(IBB aCEAcE/2(gBB)M2,CEMl,BB 

M2,BB ^-1 

Af2,CE 



(10) 



We may combine equations (6) and (10) to produce an 
expression relating the initial separation of the binary to 
that after the common envelope phase. 

acE 

Oi 



-Fi^CE = 



= -^i— *BB ' JbB^CE 



(11) 



The post-common-envelope system will consist of the 
white dwarf and the helium star in a tight binary. If the he- 
lium star is sufficiently massive, it will explode as a super- 
nova, producing a neutron star. The instantaneous mass-loss 
together with the kick the neutron receives at birth will un- 
bind some of these binaries, while others will remain bound. 
By sampling the neutron star kick distribution of Hansen 
& Phinney (1997), and allowing for the effects of mass-loss 
during the supernova, we are able to produce a population 
of WD-NS binaries and plot their binary properties and 
system kick velocities, given a population of pre-supernova 
binaries. 

Before we inspect the population of WD-NS systems 
produced in the scenario outlined above, we must con- 
sider the constraints on the properties of the original pre- 
evolution binaries. 



2.1 Constraints on initial masses 

We require that the mass of the primary should be insuf- 
ficient to produce a supernova on its own. In other words 
the mass of the compact object produced after the case BB 
mass transfer (as given in equation [4]) does not exceed the 
Chandrasekhar mass: 



Mwd = Mco(aMy < M CH 



(12) 



By definition the primary mass exceeds that of the sec- 
ondary, ie 



Mi i > M 2 ,i 



(13) 



In order to avoid the onset of a delayed dynamical instability 
when mass is transferred from the primary to the secondary 
we require 



M 2 ,i £ Mi,i/3 



(14) 



In addition, we require that the secondary does ultimately 
produce a supernova. In other words, Mco(M2 i ce) > Mch- 



As M 2 ,ce = aMj iBB = a(Mi,i + M 2 ,i 
ment yields 

/ M"He(Mc H ; 



Mwd) 6 , rearrange- 



M 2 ,i > 



l/b 



+ Mwd - Mi,i 



(15) 



The constraints on the initial masses given in equations 
(12) - (15) are illustrated in Fig. 1. The maximum allowed 
mass for the primary is set by equation (12) and is given by 
Mi,i = 8.3Mq. The upper mass limit for M2,i is given by 




Figure 1. Plot of initial primary mass Mi as a function of sec- 
ondary mass M 2 (both in solar units) showing constraints for the 
formation of a WD-NS binary as described in section 3. 



equation (13), while the lower-mass limit is given by equa- 
tion (15) for lower values of Mi,i and by equation (14) for 
larger values of Mi,i. 

2.2 Constraints on initial separations 

The evolutionary pathway described above places a number 
of constraints on the initial separation, di, of the binary if a 
WD-NS system is to be produced. 

For the first phase of mass transfer to be radiative Case 
B mass transfer, we require that the radius of the primary 
when it first fills its Roche lobe, Ri,i, is larger than the 
minimum radius required for Case B transfer, i? m in,B(Mi i i), 
and smaller than the maximum radius allowed for radiative 
case B mass transfer, i? m ax,rB(Mi j i). In other words 



(16) 



flmin,B(Ml,i) fimax,rB(Ml,j) 

/i(«0 HqO 

where Ri ; i = fi(qi)cn with q\ — Mi,i/M2,i and fi(q) = 
0.49<T 2/3 /(0.6<T 2/3 + ln(l + g~ 1/3 )). In order for Case BB 
mass transfer from the primary to occur, we require the he- 
lium star to fill its Roche lobe. In other words 

o, < fi r ,H ; (Ml \ b) ( i7 ) 

Fi_B/l(?B) 

where qs = M13/M23. In order for the late (convective) 
Case B mass transfer from the secondary to occur, leading 
to the onset of a common envelope phase, we require that 
the secondary fill its Roche lobe, but not until it has a radius 
larger than the maximum allowed for early (radiative) Case 
B mass transfer, ie 



Rn 



3(M 2 ,bb) 



^*i— >Bb/2 (QBb) 



< Oi < 



-Rmax,B(M2,BB) 
^i— »Bb/2 (QBB ) 



(18) 



where q B B = Mi iB b/M2,bb and -R m ax,B(M 2 ,BB) is the max- 
imum radius allowed for Case B mass transfer for a star of 
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mass M2,bb- A CE envelope can also follow from a donor 
with a radiative envelope if the mass ratio is big enough for 
a delayed-dynamical instability to occur. The systems of in- 
terest here have a white-dwarf primary mass of 1—1.4 Mq 
and a donor star whose mass is above the lower limit for 
evolving to a supernova, ie M2 > 8.3 Mq. This means that 
the lower limit implied by equation (18) will not be a restric- 
tion to the onset of a CE phase and thus on the formation 
of WD-NS binaries. 

We require that the helium main-sequence-star sec- 
ondary does not fill its Roche lobe at the end of the common 
envelope phase. 



(to > 



-RhoMs(M2,Ce) 

-F1^ce/2(<7ce) 



(19) 



where ifaeMS (M2,ce) is the radius of the helium main- 
sequence-star secondary and §ce = Mi.ce/A^ce- 



3 EVOLUTION IGNORING THE HE STAR 
RADIUS 

We now consider the evolution of binaries following the sce- 
nario descibed above, applying the constraints given in equa- 
tions (12) - (19). The evolution of the binary separations and 
stellar masses were computed using equations (1) — (11). The 
functions i? min ,B(M), i? ma x,rB(M), and i? max ,B(M) were 
tabulated from Bressan et al. (1993), Rmnx,He(M) from 
Paczyhski (1971) and Habets (1986), and M CO (M) from Ha- 
bets (1986). Linear interpolation was used to evaluate them 
for a particular mass, M. We note that the only informa- 
tion which we draw from Paczyhski (1971) is i? m ax,Ho(M), 
and this only to compare with results of Habets (1986), on 
which we mostly rely. The radius of a zero-age helium main- 
sequence star, RhcMs(M), can be approximated as follows 
(Habets 1986) 



log 



-RhoMS 

R<7) 



-0.68 + 0.67 log 



(20) 



For illustration we consider here systems having an ini- 
tial primary mass, Mi,i = 7.25Mq. We follow the evolu- 
tion of binaries having a range of initial separations, a ; , and 
secondary masses, Mi t \ (within the range allowed by equa- 
tions [13] - [15]). We thus synthesize a population of bi- 
naries containing a white dwarf and a helium star on the 
verge of exploding as a supernova. We then apply a ran- 
dom kick to the neutron star formed in the supernova explo- 
sion (drawn from the kick distribution of Hansen & Phinney 
[1997]) and combine it with the dynamical consequences of 
the instantaneous ejection of the helium star envelope. A 
fraction of the binaries will be broken up by the supernova; 
those that remain form our population of WD-NS systems, 
and should contain systems resembling both J1141-6554 and 
B2023+46. 

The initial parameter space (ie Oi and M24) which pro- 
duces WD-NS systems is shown as a dotted region in Fig. 
2, where the five lines drawn are the constraints given in 
equations (16) - (19). Here we have assumed qceAce = 0.5. 
The range of values of M2 is constrained by equations (13) 
and (14). The top line in Fig. 2 is the upper limit given in 
equation (17) (labelled U17) and the bottom line is the the 
lower limit from equation (16) (labelled Li6) with other lines 




Figure 2. Plot of initial separation a; as a function of secondary 
mass M2(both in solar units) showing constraints for the forma- 
tion of a WD-NS binary as described in section 3. 
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Figure 3. Plot of eccentricity e as a function of separation, a psr 
(in solar radii) for WD-NS systems produced through the evolu- 
tionary scenario described in section 3. The two open circles are 
the observed systems J1141-6545 and B2023+46. 



being labelled accordingly. We therefore see from the figure 
that the lower limit on initial separations is constrained by 
the lower limits of equations (16) and (19), ie the require- 
ments that the first phase of mass transfer is Case B, and 
that the helium main-sequence star secondary does not fill 
its Roche lobe at the end of the common envelope phase. 
The upper limit on initial separations is given by the upper 
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limit from equation (16) for M 2 < 4 Mq, ie by requiring the 
first phase of mass transfer from the primary to occur while 
the star is still radiative. For more massive secondaries, the 
upper limit on initial separations is given by the upper limit 
from equation (18), by requiring that the secondary fill its 
Roche lobe causing the onset of a common envelope phase. 
The lower limit from equation (18) divides systems enter- 
ing a common envelope phase via convective mass transfer 
(above Lis) and via a delayed dynamical instability (below 
Lis). 

The resulting population of WD-NS systems is shown 
in Fig. 3, where we plot the semi-major axes of the binaries, 
ttpsr, as a function of their eccentricities, e. The observed 
locations of J1141-6545 and B2303+46 are also shown as 
open circles (separations of 1.89 R© and 31.0 R© respec- 
tively, where we have assumed total masses of 2.3 M© and 
2.64 M© respectively). We see from this plot that under 
the assumptions of section 2, we can easily produce J1141- 
6545, but that B2023+46 is wider (by a factor of about 
2) than any system of a similar eccentricity produced by 
our population synthesis code. Wider systems will be pro- 
duced if we increase the value of oceAce to values in excess 
of unity, but then this makes it much harder to produce 
J1141-6545. In fact, in this approach, one cannot have both 
systems explained with the same value of qceAce- Whether 
acEAcE <i 1 implies o?ce > 1 depends on the value of Ace, 
which is uncertain, although values greater than unity seem 
to be possible (see Dewi & Tauris 2000, and Tauris & Dewi 
2001). 



4 ROLE OF THE HELIUM STAR RADIUS 
4.1 B2023+46-like evolution 

We consider now the evolution of the helium main-sequence 
star after the common envelope phase. In the previous sec- 
tion we assumed that its mass and radius did not evolve 
prior to the supernova explosion. Neither assumption is cor- 
rect. In isolation, the helium star will evolve to a radius 
i?max,Ho(M2,cE). If this star is to avoid filling its Roche lobe 
and transferring material to the white dwarf, we require 



cii > 



R„ 



,He(M 2 , C 



Fi^CE/2(<?CE) 



(21) 



where qcE = -M^ce/A^ce- The initial parameter space (ie 
(Zi and M2,i) which satisfies equations (16) - (19) and equa- 
tion (21) is shown in the dotted region in Fig. 4, where we 
have assumed o:ceAce = 0.5. Note that the initial phase 
space is much reduced compared to Fig. 2. If equation (21) 
is satisfied, the helium star will lose its envelope as a wind 
before the remaining core explodes as a supernova, where 
the pre-supernova core mass is given by 



Mi- 



= M C o(M 2 , 



(22) 



This mass loss will also affect the separation of the two stars. 
The relative change in separation will be given by 



-FcE-»PSN 



apsN Mwd + Af 2j cE 



acE Mwd + -M2 ,psn 



(23) 



Hence the mass loss will increase the separation of the 
two stars prior to the supernova. Continuing the evolu- 
tion of the binaries as before we produce a population of 




0.6 
log(M 2 ) 

Figure 4. Plot of initial separation a; as a function of secondary 
mass M 2 (both in solar units) showing constraints for the forma- 
tion of a WD-NS binary as described in section 4, where the 
actual radius of the pre-supernova helium star is considered. 



post-supernova systems, as illustrated in Fig. 5. Note from 
this figure that B2023+46-like are on the edge (but within) 
the phase space of allowed systems. No systems similar to 
Jl 141-6545 are produced. It is in fact quite impossible to 
produce a system resembling J1141-6545 without the he- 
lium star filling its Roche lobe in the post-common-envelope 
phase. The timescale for a wider, more eccentric system to 
evolve into a J1141-6545-like system via the emission of 
gravitational waves is prohibitively long. We therefore con- 
clude that this type of evolution is able to explain the pro- 
duction of B2023+46 but is unable to explain the origin of 
J1141-6545. 



4.2 J1141-6545-like evolution 

We consider next the case where the helium star fills its 
Roche lobe prior to the supernova, ie 



i?max,Ho(A/'2 I CE) 
Fi^CE/2(<7CE) 



(24) 



where qcE = Mi j ce/M2 i ce- A number of outcomes are then 
possible. If a second common envelope phase were to en- 
sue, the core of the helium star and the white dwarf would 
probably merge and we would be left with a single object 
rather than a binary. Alternatively, one might imagine a 
phase of mass transfer similar to that envisaged for Cygnus 
X-2 (King & Ritter 1999, Kolb et al. 2000) where mate- 
rial was transferred rapidly from a donor on to a compact 
object (in that case, a neutron star) but rather than being 
accreted, the vast bulk of the material was ejected from the 
system. Taking the data from Habets (1986), we estimate 
that Mkh = M/tkh for He-stars in the mass range of inter- 
est ~ 2.5M©- ~ 4M© is of the order of ~ 10 _3 M©yr _1 if 
the radius of the He-star is Rh c ~ 3i?© .At the same time the 
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1000 



Figure 5. Plot of eccentricity e as a function of separation, a pBT 
(in solar radii) for systems produced without mass transfer after 
the common-envelope phase as described in section 4. As in Figure 
3, the two open circles are the observed systems J1141-6545 and 
B2023+46. 




Figure 6. Plot of eccentricity e as a function of separation (in 
solar radii) for WD-NS systems undergoing mass transfer after 
the common envelope phase as described in section 4. As in Figure 
3, the two open circles arc the observed systems J1141-6545 and 
B2023+46. 



Eddington accretion rate for a WD of mass Mwd = 1.2Mq 
is M E dd ~ 6.4 l(n 6 M yr _1 ; therefore, M KH /M E dd ~ 10 2 
for the situation of interest. Since mass transfer starts from 
the more massive component, the actual mass transfer rates 
are likely to peak at even higher rates. Assuming the binary 
follows this Cygnus X-2 like evolution, and that the ejected 
material carries with it the specific angular momentum of 
the white dwarf, we have the following expression for the 
relative change in separation due to mass transfer from the 
Roche-lobe-filling helium star on to the white dwarf 



i*CE-»PSN = 



M 2 ,ce \ 2 Mwd + M 2 , CE s 
M 2j psn / Mwd + M 2j psn 



exp 



2(M 2 ,psn - M 2 , C e) 
Mwd 



(25) 



where M 2j psn is evaluated using equation (22). After this 
phase of mass transfer the core of the helium star explodes 
as a supernova. Again we synthesize a population of post- 
supernova binaries which is shown in Fig. 6. It is clear from 
this figure that J1141-6545-like systems can be produced 
via this route but that B2023+46-like systems cannot be 
a product of systems which have undergone mass transfer 
after the common envelope phase. 

We therefore see two distinct evolutionary paths: sys- 
tems undergoing mass transfer after a common envelope 
phase produce systems like J 114 1-6545 whilst those like 
B2023+46 are produced when the pre-supernova helium 
star fails to fill its Roche lobe. All computations of He 
star evolution agree on the fact that for single He stars 
i?max,He(M) drops rapidly from w 10 2 R to a few Rq be- 
tween Mho ~ 2.5 M and w 4M . It is this rapid drop in 
i H (M) which leads to the clear division between the 
evolutionary channels for J1141-type and B2303-type sys- 
tems. Where exactly this dividing line is is less important 
than that it exists. 



5 DISCUSSION 

5.1 System birthrates 

We now consider the production of J1141-6545-like and 
B2023+46-like objects from binaries containing a range of 
primary masses. As can be seen from Fig. 1, constraints on 
the initial masses of the two stars (as decribed in section 
2.1) restrict the primary mass to 4.6M ^ Mi ^ 8.3M . In 
Fig. 7, we show the computed fraction of binaries produc- 
ing J1141-6545-like and B2023+46-like objects, assuming 
the binaries to be distributed uniformly in log(ai/R ) and 
1 < l°g( a i) < 3. From this figure we note that the further 
restrictions on the parameter space produce the two types of 
binary (as described in section 2.2) are such that no binaries 
of either type are produced for Mi <J 6.4Mq. This in turn 
implies that the mass of white dwarfs generated in this way, 
Mwd i<; 1 Mq . We note also that the fraction of binaries pro- 
ducing J1141-6545-like systems is far higher than that for 
B2023+46-like systems. As also shown in Fig. 7, we consid- 
ered different distributions of secondary masses, ie different 
values of a where cW 2 /dM 2 tx M^ a . A flat distribution (ic 
a = 0) produced a larger fraction of both types of binary 
compared to a Salpeter-like distribution (ie a = 2.35). This 
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M, 

Figure 7. The fraction of binary systems, /, producing either an 
J1141-6545-likc system (solid lines) or a B2023+46-likc system 
(dashed lines). In each case the lower of the two lines assumes a 
Salpcter-like distribution of secondary masses, whilst the upper 
line assumes a flat distribution of secondary masses. 



a 




Nun 


^2303 


2.35 


8.75 


xl0~ 5 


2.45 xlO" 6 


0.00 


6.20 


xlO- 4 


5.70 xlO" 5 


-2.35 


7.60 


xlO- 4 


1.17 xlO" 4 



Table 1. The galactic production rate (in units of systems yr x ) 
of J1141-6545-like and B2023+46-like systems (JVn 4 i and 7V 2 303 
respectively) as a function of distribution of secondary masses, a, 
where dN 2 /dM 2 oc M^ a . 



is because larger secondary masses are required to produce 
either J1141-6545-like or B2023+46-like systems. 

We can now compute the galactic formation rate of sys- 
tems assuming a primary formation rate dN\ = M 1 _2 3B dMi 
(which is equivalent to a galactic star formation rate of 
~ 3M0/yr). We also assume that one quarter of all mas- 
sive stars are formed in binaries having the separation range 
1 < l°g( a < 3. Hence the formation rate of systems of 
type X, derived from binaries having an original primary 
mass Mi, is dNi t x = fxdNi/4, where fx is the fraction 
of binaries, having the separation range 1 < log(ai) < 3, 
which produce systems of type X (as shown in Fig. 7). In 
Table 1 we list the galactic formation rate of systems as a 
function of a, for a co A = 0.5. We note that for both J1141- 
6545-like and B2023+46-like systems, the formation rate 
increases with decreasing a, as a larger fraction of the origi- 
nal binaries contain more massive secondaries. The relative 
formation rates of the two types of systems is also a function 
of a as B2023+46-like systems require a secondary of mass 



0.001 u i i i 




Figure 8. The galactic production rate (in units of systems yr x ) 
of J1141-6545-likc systems (solid lines) and B2023+46-like sys- 
tems (dashed lines) as a function of ct co A. In each case the lower 
of the two lines assumes a Salpeter-like distribution of secondary 
masses, whilst the upper line assumes a flat distribution of sec- 
ondary masses. 



similar to that of the primary (as shown in Fig. 4) whereas 
J1141-6545-like systems are produced for a broader range 
of secondary masses. We now consider the effect of changing 
the value of a cc A. In Fig. 8 we plot the galactic produc- 
tion rates for J1141-6545-like and B2023+46-like systems, 
for the two cases a = 2.35 (ie Salpeter-like distribution of 
secondary masses) and a = 0.0 (ie a flat distribution of 
secondary masses). The production rate of J1141-6545-like 
systems is relatively independent of a cc A whereas the pro- 
duction rate of B2023+46-like systems increases with a co A. 
The formation rate of J1141-6545-like systems appears to 
be larger than for B2023+46-like systems for all reasonable 
values of a and q co A. 

We now consider the formation rates inferred from the 
two observed systems. J1141-6545 is at a distance ~ 3 kpc 
and has a characteristic age derived from the pulsar spin- 
down of 1.4 Myr. The distance to B2023+46 is less well- 
known but is thought to lie between 3-10 kpc. It has a 
spin-down age of 30 Myr. The beaming of pulsar emission 
means that we will only see those which are pointing towards 
us. We might only see 1/5 of the actual population. Including 
the effect of beaming, the distances to both systems would 
suggest that there are currently 50 - 500 similar systems 
thoughout the Galaxy. From the pulsar characteristic ages, 
this in turn implies formation rates of 5 x 10~ 5 — 5 x 10~ 4 
yr" 1 for J1141-6545-like systems and 2 x 10~ 6 - 2 x 10~ 5 
yr -1 for B2023+46-like systems. These values are consistent 
with the range seen in Table 1. From Figure 8, we note that 
in addition that the observed formation rates of B2023+46- 
like systems would seem to suggest that q co A <; 0.3. We 
therefore conclude that the scenario described in this paper 
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Figure 9. The kick velocity received by the WD-NS binaries Vbj n 
(in km/s) as a function of their semi-major axes a psr (in solar 
radii). Crosses are for systems produced via the J1141-6545-like 
route, whilst the open squares are for systems produced via the 
B2023+46-likc route. 

is able to account for the two observed systems PSR J1141- 
6545 and B2023+46. 

5.2 Kick velocities 

We consider now the system kick received by J 114 1-6545- 
like and B2023+46-like binaries. The system kick of each 
binary is a combination of the kick the neutron star receives 
at birth with the kick resulting from the mass-loss in the 
supernova explosion of the helium star. The distribution of 
system kick velocities are shown in Fig. 9. From this fig- 
ure we note that the kicks span a large range of values, the 
largest being likely to remove the system from the Galaxy 
entirely (when combined appropriately with the initial or- 
bital velocity of the pre-supernova binary). In a majority of 
cases, such kicks will at least take the binary significantly 
away from the plane of the galactic disk (where the sys- 
tems will have been formed). The kicks for systems most 
closely resembling J1141-6545 range between 50 km/s and 
300 km/s. For B2023+46 the equivalent range is 100 km/s 
to 200 km/s. These are both consistent with the observed 
positions of J1141-6545 and B2023+46. J1141-6545 has a 
galactic latitude of —3.86 degrees which, at a distance of 3.2 
kpc and assuming an age of 1.4 Myr, implies a velocity out 
of the galactic plane <J 150 km/s. B2023+46 has a galactic 
latitude of —12.02 degrees which, at a distance of 3 - 10 kpc 
and assuming an age of 30 Myr, implies a velocity out of the 
galactic plane of ~ 20 — 60 km/s. 

5.3 Merger rates 

In this section we discuss the subsequent evolution of the 
J1141-6545-like objects. As can be seen from Fig. 6, in these 



Figure 10. The cumulative distribution of timescales for the pop- 
ulation of inspiralling WD-NS binaries to come in to contact via 
the emission of gravitational radiation, t CO ntact (in years). 

systems, the white dwarf and neutron star have separations 
~ 1— 1OR.0, and a wide range of eccentricities. These systems 
will spiral in as they lose angular momentum and energy via 
the emission of gravitational radiation. One may compute 
the subsequent evolution of the binary separation, a, and 
eccentricity, e, using the following expressions (Peters 1964) 

da 64G ,3 MwdM ns (Mwd + M NS ) 

dt ~ 5c 5 a 3 (l-e 2 ) 7 /2 X 




de _ 304G 3 MwdM ns (Mwd + M NS )e / 121 2 \ 

dt ~ 15c 5 a 4 (l - e 2)V2 x V 1 + 3 04 e J [ ' 

We thus computed the time required for the binaries to 
evolve into contact. The cumulative distribution for contact 
times is shown in Fig. 10. We see from this figure that over 
half of the systems will merge in a timecale Ss 10 8 years, 
and more than 95% within a Hubble time. These merging 
systems may make an important contribution to the popu- 
lation of gamma-ray burst progenitors as will be explored in 
a later paper. 



6 SUMMARY 

We have described an evolutionary pathway for producing 
white dwarf-neutron star binaries where the white dwarf is 
made first. This explains the observed systems J 114 1-6545 
and B2023+46. In this scenario, the original primary trans- 
fers its envelope to the secondary conservatively during a 
radiative Case B phase of mass transfer leaving only the 
helium-star core. This helium star then evolves, filling its 
Roche lobe and leading to the second phase of mass transfer 
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(Case BB) where the envelope of the helium star is trans- 
ferred to the secondary (which is still on the main sequence) . 
In the process the primary becomes a CO or ONeMg white 
dwarf. The secondary evolves into contact on its nuclear evo- 
lutionary timescale. If the system is sufficiently wide, this 
third phase of mass transfer leads to a common-envelope 
phase during which the white dwarf and the helium-star 
core of the secondary spiral together as the common enve- 
lope of gas is ejected from the system. The post-common- 
envelope system consists of the white dwarf and the helium 
star in a tight binary. Providing sufficient mass has been 
transferred to the secondary during earlier phases in the bi- 
nary evolution, the helium star may be massive enough to 
produce a neutron star via a supernova explosion. The post- 
common-envelope but pre-supernova evolution depends on 
the separation of the binary. 

Systems like B2023+46 are produced in relatively wide 
systems, where the helium star avoids filling its Roche lobe 
as it expands. Rather it loses mass via a wind. The core of 
the helium star eventually explodes as a supernova produc- 
ing a neutron star. Assuming the neutron star receives a kick 
at birth, the majority of these systems are broken up, pro- 
ducing high-velocity, single, white dwarfs and neutron stars. 
In a minority of cases, the neutron star and white dwarf 
produce a binary having a separation ~ 5 — 50R© . 

For closer binaries, the helium star in the post-common- 
envelope system fills its Roche lobe. In this paper, we have 
proposed that a phase of mass transfer follows where the 
helium star envelope is transferred to the white dwarf at 
highly super-Eddington rates. Rather than being accreted, 
we propose that this material is ejected from the system. 
Subsequently the helium-star core explodes as a supernova, 
producing a neutron star. As these binaries are significantly 
tighter than those producing B2023+46-like systems, the 
vast majority remain bound, having separations in the range 
~ 0.5 — 5Rq. This is the way we believe J1141-6545 was 
formed. 

We have computed the expected formation rates for 
J1141-6545-like and B2023+46-like systems and shown 
that they are both large and consistent with the two ob- 
served systems. 

We have shown that the vast majority of J1141-6545- 
like binaries will spiral in and merge in less than a Hubble 
time due to the emission of gravitational radiation. Given 
their relatively high formation rate (~ 1CP 4 — 1CF 5 yr _1 ), 
these systems may make an important contribution to the 
population of gamma-ray burst progenitors. 
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